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We present a metamaterial-based random polarization control plate to produce incoherent laser 
irradiation by exploiting tiie ability of metamaterial in local polarization manipulation of beam 
upon transmission via tuning its local geometry. As a proof-of-principle, we exemplify this idea nu- 
merically in a simple optical system using a typical L-shaped plasmonic metamaterial with locally 
varying geometry, from which the desired polarization distribution can be obtained. The calcu- 
lating results illustrate that this scheme can effectively suppress the speckle contrast and increase 
irradiation uniformity, which has potential to satisfy the increasing requirements for incoherent laser 
irradiation. 

PACS numbers; 



I. INTRODUCTION 

Coherence is a fundamental property of laser, which 
facilitate laser to be widely applied in those fields re- 
quiring high coherent light sources. However, just as a 
coin has two sides, coherence is harmful to some laser 
application fields. High coherence usually results in un- 
wanted speckle noise in laser display (or laser projection 
imaging or laser TV), thereby decreasing the imaging 
resolution [1]. In laser fusion and laser heat processing, 
coherence makes the intensity distribution on the focal 
plane not uniform enough In some spectroscopy ex- 
periments, high coherence laser is also undesirable. Thus, 
in these applications, it is urgently need to eliminating 
the laser coherence. 

In essence, controlling the polarization or phase of laser 
beam could achieve incoherent laser irradiation. The key 
points of beam docoherence through polarization control 
can be summarized as follows. It first breaks a laser beam 
into two or many beamlets with orthogonal or random 
polarization distribution and then converge them; the ir- 
radiation on the focal plane is thereby incoherent. Some 
concrete schemes have been proposed, such as using a 
birefringent wedge to create two orthogonally polarized 
beams with a selected angular separation and employing 
an optic made by wave plates to scramble the polariza- 
tion distribution in the near field Generally, these 
polarization control devices are made of nematic liquid 
crystal, KDP crystals, and crystal half-wave plates, etc. 

On the other hand, as artificial periodic subwavelength 
structures, metamaterial can exhibit strong anisotropy 
by suitably designing its microscopic structure unit 0,[3- 
Generally, naturally existing anisotropic mediums exhibit 
very small birefringent index in optical domain, which 
need relatively thick slabs to achieve a polarization con- 
trol device (e.g., hundreds or thousands of wavelengths 
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for a half- wave plate) . However, it is possible for metama- 
terial to achieve a polarizer with the same polarization 
conversion efficient as convention anisotropic mediums, 
with its thickness less than a wavelength and transverse 
dimension in the order of wavelength |9l-[l3j . So it holds 
great potential in the future nanophotonics applications, 
and is particularly amenable to miniaturization. From 
this point of view, it would be interesting to apply meta- 
material in polarization control devices. 

In this work, a metamaterial-based random polariza- 
tion control plate (MRPCP) is proposed for eliminating 
the laser coherence and producing incoherent laser irra- 
diation. Locally tailoring the geometric structure param- 
eters of an anisotropic metamaterial, we could obtain a 
locally varying polarization response of laser beam pass- 
ing through it, thereby scrambling the polarization state 
of beam in the near field and finally eliminating the laser 
coherence on the focal plane. Without loss of generality, 
we will exemplify our scheme numerically using a varia- 
tion of the typical L-shaped metamaterial, although the 
approach can be applied to an arbitrary metamaterial 
geometry offering enough number of free design parame- 
ters. The calculating results show that the MRPCP can 
effectively suppress the speckle contrast and increase ir- 
radiation uniformity. 



II. THEORY AND SCHEME 

The mechanism of MRPCP can be considered by em- 
ploying the multi-beam interference formula 

/(r,t) = 

i 

-f ^ 2 [a, ■ aj] E,Ej cos[$,(r, t) - $j (r, i)](l) 

where Oi^j, Ei^j, and represent the polarization direc- 
tion vector, field amplitude, and phase of the (z, j )th field, 
respectively, /(r, t) is the total intensity with location (r) 
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FIG. 1: Schematic picture of the MRPCP. Each square repre- 
sents a elements. The arrows indicate e-rays of the anisotropic 
metamaterial elements. 



and time (t) dependence. It is obvious that incoherent 
addition of the beamlets is achieved when their polariza- 
tion vectors are random, since the vector summation of 
the interference term (the second term) in Eq. ([T]) is zero. 
This is the basis of the polarization docoherence. 

The intensity distribution of MRPCP combined with 
a rando m p hase plate (RPP) is referred to as laser 
speckle [l4 [lH . The statistical property of the speckle 
intensity can be considered as the addition of two or- 
thogonally polarized speckle patterns (such as in x- and 
y-directions) [H, [13] ■ Due to the non-equivalence of the 
two speckle intensities {Ix 7^ ly), the probability density 
of the intensity can be calculated by [3) US] 
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with Ix^y indicating the average intensities of the two 
speckle patterns. Because of the incoherency of the two 
orthogonal components, we have Itotai = Ix + ly, and 
then Eq. ([2]) can be written as 
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where Ai and A2 represent Ix/ Itotai and ly/ Itotai, re- 
spectively. The speckle contrast is also reduced and de- 
termined by the ratio of Aiand A2. The formula for 
calculating the contrast is jl7| 



^Al + Al/{A, + A2). 



(4) 



When Ai = A2, C take the minimum value 0.707. It 
means that the speckle contrast can be reduced 29.3% at 
most by using MRPCP. 

For the convenience of calculation, we now design a 
MRPCP using L-shaped plasmonic metamaterial, which 



requires it producing a random polarization distribution. 
A simple scheme is to make the optic-axis of the meta- 
material taking random orientations for different local 
regions, namely, locally varying the structure geometry 
could create an MRPCP. The schematic picture of the 
MRPCP is shown in Fig.[l]where the arrows indicate the 
e-rays of the anisotropic metamaterial elements. It is as- 
sumed that the number of the elements is to x n (m and 
n are both positive integers). An element is comprised 
of several L-shaped unit cells with the same orientation 
(e-rays) and subwavelength dimension, and the number 
of the unit cell for each element is not restricted. 

We first consider the polarization change of an element. 
The schematic picture of the metamaterial unit cell is 
shown in Fig. [21(a). The "L" is assumed to be made of 
gold with 20 nm thickness and built on a dielectric sub- 
strate {us = 1.45). They together form a unit cell of the 
metamaterial. Other structure parameters are labeled in 
the figure. We start by employing the complex Jones ma- 
trix of this metamaterial element [31 ■ It can be utilized 
to describe the transmission of coherent light through any 
dispersive optical system. The complex Jones matrix T 
associating the complex amplitudes of the incident field 
(Eix.y) with the transmitted field {Etx,y) is given as 
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Here, it is assumed that 6 — and the incident plane 
wave propagates in positive z-direction (normal inci- 
dence) . 

Both arms of the "L" are mutually perpendicular 
with identical structure dimensions. This indicates that 



Txx — 



I-'yy and Txy 



The arms can be seen as 



oscillators with respective eigenfrequency, which are as- 
sociated with the excitation of carriers representing the 
free-electron gas of the metal. The carrier dynamics of 
one arm is affected by the external electric field as well as 
the conductive coupling from another arm [l^ [23| . Due 
to the identical structure dimensions, the eigen resonant 
responses of the two arms are identical as well as the con- 
ductive coupling responses. So we can safely obtain that 
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FIG. 2: (a) Schematic of the unit cell with the geometrical 
parameters: b = 250 nm, c = 60 nm, and the thickness of the 
substrate d = 100 nm. The lattice constant along both x- and 
y-directions are / = 500 nm. The thickness of the gold "L" is 
20 nm and not shown in this figure, (b) Rotating of "L" from 
the X — y coordinate system to u — v coordinate system. 
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of freedom to design the MRPCP, i.e., changing the ge- 
ometric dimensions of the L's arms. Arms with different 
length and thickness wih induce different optical reso- 
nance responses (i.e., different T matrix elements). The 
T matrix can be simulated for a number of slightly vary- 
ing dimensions of "L" , resulting in a look-up table linking 
geometric parameters with T matrix elements. This ta- 
ble represents the set of possible building blocks of the 
MRPCP. 



FIG. 3: Magnitude of the transmission coefficient versus fre- 
quency /. 



-^xx — yy "'^^^ xy — ^yx- 

In order to solve T^x and Txjyj we use the finite- 
difference time-domain method [21] to simulate the meta- 
material unit cell shown in Fig.[2Ja). In the simulations, 
the metal permittivity is described by Drude model: 
e(a;) = f — oj^/iuP' + i'jui), with the bulk plasma fre- 
quency ujp and the relaxation rate 7. For gold, from 
literature [2^, tOp = 1.37 x IQ-* THz and 7 = 40.84 THz. 
The magnitudes of T^x [Tyy) and T^y {Tyx) are plotted 
in Fig. [3]when the incident electric field only polarized in 
X- (?/-)direction. Each curve has two valleys (or peaks) 
which means that each oscillator has two distinct eigen- 
modes near the frequencies 173 THz and 290 THz, re- 
spectively [l^l. It also indicates that, if the incident elec- 
tric field is polarized in either direction of the two arms, 
cross-polarized radiation will generate and of course re- 
sult in polarization change in the transmitted field. In 
fact, the structure dimensions including the length and 
thickness of the arms can affect the resonant responses of 
the metamaterial, e.g., longer length or thinner thickness 
of the arms would decrease the resonant frequencies. 

Now, the Jones vector of the transmitted field can be 
calculated according to Eq. dS]). Without loss of gener- 
ality, when takes an arbitrary value [see Fig.[2I[b)], the 
Jones vector can be obtained by transformed the x — y 
coordinate system to the u — v coordinate system: 

(t:)^«<'"(r::T::)"'-'"(t:)- 

where 

/cose -sin^y 

According to Eq. ([5]), at any frequency, we can obtain the 
Q dependence of the polarization states of the transmitted 
field. It's worth noting that the T matrix can be diago- 
nalized by a rotation of 7r/4 or -7r/4 since T„„ — T^y and 
Tuv = Tyu, which indicates that the angle bisector direc- 
tion of "L" and its orthogonal direction represent the e- 
and o-rays of the anisotropic metamaterial, respectively. 

If both arms of the unit cell are non-identical, T^u is 
not equal to Tyy, but Tyy = Tyy can still be obtained. 
The T matrix could no longer be diagonalized whatever 
value 9 is. Therefore, actually, we have another degree 



III. APPLICATION OF THE MRPCP FOR 
SUPPRESSING THE SPECKLE CONTRAST 
AND INCREASING IRRADIATION 
UNIFORMITY 

As a proof-of-principle for the effectiveness of MRPCP 
in suppressing the speckle contrast and increasing irra- 
diation uniformity, a simple optical system is considered 
here and schematically illustrated in Fig.|4l A plane wave 
Ei{x,y) polarized in the x-direction {Eiy = 0) with unit 
amplitude is assumed for incident laser. It illuminates 
the MRPCP and is converged by a lens after passing 
through a RPP. 

According to Eq. we can obtain the transmitted 
field of X and y components from the (m, n)th metama- 
terial element of the MRPCP: 

/ Etx \ _ ( {Tyy — Tyy SlU 26) \ , , 

\Ety)^\ E,,TyyCOs20 )■ 

From Eq. ([5]), it is easily to find that light arising from 
different metamaterial elements with random 9 values 
will undergo different polarization change. When a laser 
beam passing through the random MRPCP, its trans- 
verse distribution of polarization state will be scrambled 
in the near field and thus make its spatial coherence be- 
come worse. Then, after passing through a RPP adjacent 
to the MRPCP, we converge the beamlets by a lens with 
focal length /q. As Tyy and T„„ have been calculated 
above by the finite-difference time-domain method, the 
field distribution on the focal plane can be derived from 




laser MRPCP RPP lens focal plane 

FIG. 4: Schematic picture of the simplified optical system. 
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FIG. 5: (Color online) Normalized one-dimensional intensity 
distribution (cross section at t/ = 0) of the speckle pattern for 
(a) only using RPP and (b) MRPCP+RPP. The green solid 
lines represent the intensity envelope. 



the Fresnel-diffraction integral formula: [23 



exp(ifc/o) exp 
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where x' and y' are axes parallel to the x and y axes, 
respectively, SajS is the space range of the RPP illumi- 
nated by the (m, n)th element of the MRPCP, (j)Rpp is 
the random phase of the RPP element which is equal to 
either or tt, and j G {x,y). The j component of the 
total electric field is given by 



(10) 



and then the intensity of the total electric field can be 
written as 

Itotai{x',y') = L^{x',y') + Iy{x',y') 

= 5]ii?,(x',y')P + Ei^^(^''y')Rii) 



The normalized one-dimensional intensity distribu- 
tions of the speckle patterns on the focal plane for only 
using RPP and MRPCP-^RPP are shown in Fig. El The 
dense spikes come from the interference among the beam- 
lets, and the intensity envelope (the so-called Airy spot 
indicating by green lines in Fig. [5]) is determined by the 
diffraction of each element of RPP. One can notice that 
there is almost no zero-intensity in the central main re- 
gion for MRPCP-I-RPP, while numerous zero-intensities 
for only using RPP. Also shown in Fig. O the introduc- 
tion of MRPCP substantially reduces the contrast be- 
tween the "spikes" and envelope. For further understand- 
ing these results, we will statistically analyze the inten- 
sity distributions obtained using them. The statistical 
results are calculated by considering thousands of points 
of the speckle pattern on the focal plane. 

The probability density of intensity for one speckle 
pattern behaves as a negative exponential distribution 
(dashed line in Fig. 15]), such as for the case of only us- 
ing RPP [Hi. For MRPCP+RPP, statistical results from 
the focal plane show that the average intensities of the 
two orthogonal directions are non-identical [Ai k, 0.554 
and A2 ~ 0.446), because the number of the element 
of MRPCP in our real calculations is finite. Theoret- 
ically, the shape of the probability density distribution 
depends on the ratio of the two average intensities Ai 
and A2. For comparison, we directly calculate the prob- 
ability density of MRPCP+RPP from Eq. and draw 
it in Fig. [6] (blue line). It indicates that the probabil- 
ity of zero-intensity is zero and the maximum probabil- 
ity moves toward the higher intensity. This theoretical 
results are well in accordance with the statistical ones 
(indicating by green squares in Fig. |6]). The speckle con- 
trast is also reduced and determined by Eq. (jj]), here, 
C = 0.711. The reduction of the speckle contrast is about 
28.9% which is approaching to the ideal value 29.3%. 

Actually, incoherent laser irradiation are required by 
many applications fields, such as in laser display, laser 
fusion, laser heating processing, and some spectroscopy 
experiments. Our scheme holds great potential in these 
applications, especially enables to produce incoherent ir- 
radiation at very small dimensions (e.g., in the order of 



Now, we calculate the intensity distribution 
Itotai{x' ,y') on the focal plane. The parameters 
used for the calculations are a = 2 mm, /o = 1 ni. The 
operation wavelength is set as 1.5 /xm (i.e., / = 200 
THz). The total number of the metamaterial elements is 
16 X 16 whose values of 9 are random numbers between 
and 2-K produced by the computer. The number of 
the metamaterial elements should not be too few so as 
to fully randomize the polarization state. The RPP has 
6464 square phase elements (side length 0.5 mm), and is 
placed after the MRPCP. The size of the metamaterial 
elements are larger than the phase elements in order 
that one metamaterial element can illuminate at least 
one phase element. The shape and size of the speckle 
pattern are determined by those of the individual RPP 
elements. 



MRPCP-i-RPP 
(statistical) 



MRPCP-hRPP 
(theoretical) 




1 2 
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FIG. 6: (Color online) Normalized probability density 
(hotai X P{I /hotai)) of intensity distribution for only using 
RPP and MRPCP+RPP (theoretical and statistical). 
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wavelength) due to the subwavelength structure features 
of metamaterial, which is not feasible for previous meth- 
ods. An example is fabricating such structure directly at 
the tip of an optical fiber to form a compact device for 
beam decoherence [24] . 

IV. CONCLUSIONS 

We have proposed a metamaterial-based random po- 
larization control plate for eliminating the coherence of 
laser and producing incoherent laser irradiation, since 
metamaterial can offer degrees of freedom to tailor 
their transmission and polarization properties by locally 
changing its structure geometry. A random polarization 
control plate based on L-shaped plasmonic metamate- 
rial has been designed, which has locally varying struc- 



ture parameters from which the desired far-field polariza- 
tion distribution can be obtained. The calculating results 
show that it can effectively suppress the speckle contrast 
and increasing irradiation uniformity, which is potential 
to satisfy the increasing requirements for incoherent laser 
irradiation and especially enables beam decoherence at 
very small dimensions. 
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